Loop Groups and Discrete KdV Equations 



Jeremy Schiff 
Department of Mathematics, Bar-Ilan University 
Ramat Gan 52900, Israel 
email: schiff@math.biu.ac.il 

September 2002 



Abstract 

A study is presented of fully discretized lattice equations associated with the KdV hi- 
erarchy. Loop group methods give a systematic way of constructing discretizations of the 
equations in the hierarchy. The lattice KdV system of Nijhoff et al. arises from the lowest 
order discretization of the trivial, lowest order equation in the hierarchy, bt = b x . Two new 
discretizations are also given, the lowest order discretization of the first nontrivial equation 
in the hierarchy, and a "second order" discretization of bt = b x . The former, which is given 
the name full lattice KdV has the (potential) KdV equation as a standard continuum limit. 
For each discretization a Backlund transformation is given and soliton content analyzed. 
The full lattice KdV system has, like KdV itself, solitons of all speeds, whereas both other 
discretizations studied have a limited range of speeds (being discretizations of an equation 
with solutions only of a fixed speed). 



1 Introduction 



Despite the fact that numerical simulations of PDEs of KdV type can be done quickly and 
accurately these days using standard spectral methods, it is still of interest to look at discretiza- 
tions of such PDEs, and see how "integrability properties" (elastic soliton scattering, existence 
of conserved quantities etc.) are affected by discretization, and in particular to see if there are 
"integrable discretizations", that exhibit all the special properties of the underlying PDE. One 
can consider both "partial" and "full" discretizations; in the former only the spatial coordinate 
is discretized, in the latter time is discretized too. This paper focuses on full discretizations. 



1 



The difference equation usually known as discrete KdV was first studied by Hirota Hirpta. 
Using a slightly different notation from that of Hirpta, discrete KdV is the equation 

1 1 



= c(u n+ljm - u n , m+ i) (c constant) . 



l + u n+ljTn+ i l + u r 

This is a discretization of KdV, but in a rather unusual sense. The main justification for the 
name "discrete KdV" is that dh has a bilinear formulation and a family of soliton solutions very 
similar to those of KdV (see also he for rational solutions). The study of discrete KdV was taken 
further by Nijhoff and collaborators (see Nijhoff for a review and references). The work of Nijhoff 
et al. focuses on the equation 

/ 1 + fen.m+1-WA ( j + Km ~ b n+1 , m+1 \ =1 constant ) , (2) 

V p-q J \ p + q J 

which they call lattice KdV. In fact this equation is a "potential form" of discrete KdV, in the 
sense that if b rhm satisfies dp, then it is easy to check that 

bn— l,m— i b nrn 

satisfies dh with c = (p + q)/(p — q). (This is supposed to be an analog of the fact that if b(x,t) 
satisfies the "potential KdV" equation b t = \b xxx + \b 2 x + 5(t) for some function 8(t), then u = b x 
satisfies the KdV equation u t = \u xxx + 3uu x .) 

Nijhoff et aVs lattice KdV equation has an advantage over Hirota's discrete KdV in that it 
is easier to see its continuum limit (in the usual sense, to be explained shortly) as well as at 
least one nonstandard continuum limit in which it reduces to the potential KdV equation. On 
substituting p — 1/h and q — 1/k, dp becomes 

bn+l,m+l ^n,m+l ~\~ ^n+l,m ^n,m ^n+l,m+l ^n+l,m ^n,m+l ^n,m 



h k 

+ {b n ,m — &n+l,m+l)(^n,m+l — &n+l,m) = . (4) 

Taking the standard continuum limit will be taken to mean replacing 6 n m by b(x,t), b n+ i^ m by 
b(x + h,t), 6„, m +i by b(x, t + k), b n+ i jin+ i by b(x + h,t + k), expanding in powers of h and k and 
ignoring all but leading order terms. It is clear that in this limit the first term in dp2 gives —2b x , 
the second 2b t and the third 0. Thus in the standard continuum limit, lattice KdV is simply a 
discretization of b t — b x . A nonstandard continuum limit of dp2 that gives the potential KdV 
equation is as follows: Make the same replacements as before, expand in powers of h and k, but 
keep not only the leading order terms but also all terms of order h and h 2 . This gives 

( h h 2 \ ( h 2 \ 

-216, + -b xx + —b xxx \ + b t + lb t + hb tx + —b txx \ + h 2 b 2 x = (5) 
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Now write b = b — ^b x . Ignoring terms of order h 3 and above, the last equation can be written 

h = K - — i^-b xxx + -blj . (6) 

This is a "linear combination" of the flow obtained in the standard continuum limit with the 
potential KdV flow. 

The foregoing discussions raise a variety of questions. The relationship of KdV/potential KdV 
and discrete KdV/lattice KdV as it stands is rather cryptic and requires some clarification. It 
would also be good to have another integrable lattice equation from which KdV/potential KdV 
can be obtained by taking a standard continuum limit. If this is possible, then it would be good 
to know just what freedoms there are in constructing integrable discretizations. Finally, though 
this is a question that will not be addressed in the current paper, given an integrable lattice 
equation, just how much freedom is there in taking the continuum limit? 

This paper discusses the subject of discretizations of KdV using loop group methods. The 
basic fact behind the loop group approach to KdV is that the KdV equation (or, more precisely, 
the Lax pair for the KdV equation) is simply a "disguised" version of the Frobenius-integrable 
pair of linear first-order constant-coefficient ODEs 

l)u, 8 t V = l(] \)U (7) 

(here U is a 2 x 2 matrix function of x,t,l). The relation of the above system with KdV will 
be explained fully in section 2 below. In greater generality, the iVth flow (N = 1, 3, 5, . . .) in the 
KdV hierarchy is associated with the system 

W-Cl l)u, w-fi l o )"u, (8) 

which reduces to the standard system linl when N = 3. The approach proposed in this paper for 
constructing integrable discretizations of KdV is simply to discretize the system linl or lin2 (any 
explicit scheme for numerical integration of ODEs can be used) and then to apply the necessary 
"disguise" to translate this system into a discrete KdV. Section 3 is devoted to the simplest 
discretization of lin2 with N — 1, namely 

1 



+l,m 



/0 1 

I + h[ , 

V i o 



U nm) XJ ; 



n,m+l 



I + k 



I 



U n ^ m . (9) 



This is just a first-order Euler scheme with different step sizes in the x and t directions. This 
scheme gives rise to the lattice KdV equation, which, as shown above, is a first-order discretization 
of the N = 1 flow in the potential KdV hierarchy, b t = b x . As an application of the loop 
group formulation, a Backlund transformation for dp is given, and soliton solutions are derived 
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(c.f. Hirpta). A brief analysis of the soliton solutions is given, which helps clarify the rather 
schizophrenic nature of the lattice KdV equation, which on one hand is a (nonlinear) discretization 
of b t = b X) and on the other displays features of potential KdV. 



Section 4 is devoted to the simplest discretization of linl, namely 



n+l,m 



I + h 



1 

1 



u, 



n,m+l 



I + kl 







1 



U„ 



(10) 



As expected, this gives rise to a system which is, in a natural way, a first-order discretization 
of the potential KdV equation. The system is a little complicated, involving two auxiliary fields 
(reminiscent of the discretization of the sinh-Gordon equation given in bobenko), but it seems 
this is the price that has to be paid to have an integrable lattice equation that has potential KdV 
as a natural continuum limit. The Backlund transformation and soliton solutions are derived for 
this system too. 

Section 5 considers another discretization of lin2 for N — 1, namely 



+l,m 



u, 



n,m+l 



I + h 



I + k 



'0 




r 
o, 



h 2 

+ y 







1 



Un,r, 



1 



2 



hi 



h 

1 + f . 



u, 



n,m 



(11) 

(12) 



This example is worked out mainly to illustrate that the method can be extended to arbitrary order 
discretizations of linl and iin2, establishing that there is quite a lot of freedom in constructing 
integrable discretizations. Section 6 contains some concluding remarks. 

I conclude the introduction with a brief mention of some relevant literature. The approach to 
discretization taken in this paper is closely related to the approach of discretizing the scattering 
problem, which was first proposed by Ablowitz and Ladik al, and recently has been revisited by 
Boiti et al. BPPS. Several potentially interesting applications of discretizations of equations of 
KdV type have emerged recently. Nijhoff et al. Nijhoff were the first to notice the link between 
lattice KdV and the discrete conformal map equation 

(-2-n,m Zn+l,m) i,^n,m+l ^n+l,m+l) 



= S 



(s constant) , 



(13) 



\Zn,m Zn,m+l) \^n+l,m ^n+l,m+l) 

which in the case s = — 1 is a natural discretization of the Cauchy-Riemann conditions. Techniques 
related to those of this paper have been applied to dcm in dcmpaper. dcm may well play a 
significant role in the field of numerical conformal mapping. Discretizations of KdV and related 
equations have also been shown to have a role in a variety of other numerical algorithms numapps. 
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2 KdV as a linear constant-coefficient flow 



This section contains a summary of results from mel, relating the (potential) KdV equation with 
the linear constant coefficient flow linl. A rather more mathematical description can be found in 
sw. 

The general solution of linl is 

U(x,t,l)=exp(x^ i q)+^(; ljju(0,0,l) ■ (14) 

Assume that the function U(0, 0, /) is defined for \l\ = 1, and has nonzero determinant; in other 
words it is an element of the loop group LGL 2 {C) PJ3. Then evidently so is U(x,t,l). Now a 
typical element g(l) of the loop group LGL 2 (C) can be written as a product S^ 1 (l)Y(l) where 
Y{1) is holomorphic for \l\ < 1 and S(l) is holomorphic for \l\ > 1 with S'(oo) = /. This is the 
so-called Birkhoff factorization theorem, see PJ3, chapter 8. So let us write 

U(x, t, I) = S~\x, t, l)Y(x, t, I) , (15) 

(with Y holomorphic in \l\ < 1, S holomorphic in |/| > 1 and S(x, t, oo) = /) and let us try to 
find differential equations satisfied by the two "components" Y and S of U. Substituting fact 
into linl, mulitplying on the left by S and on the right by Y^ 1 gives 

-SxS^ + YsY-^S^ l^jS' 1 , -StS^ + YtY-^lS^ ^ S' 1 . (16) 

Now, if 

l\ Cl (x,t) d^t)) l 2 \c 2 (x,t) d 2 (x,t)J 
then a brief calculation shows 

/0 1\ , / -b 1\ 



\l 0J \l-v b) v 

/0 1\ / -bl-B l + v-b 2 \ t >. 

is [ Is- 1 = f )+o(r 1 ) , (19) 

\l OJ \l 2 -vl-V bl + B J V 1 

where b = —b±, v — a± — d±, B — c\ — b 2 + ai&i, V = a 2 — d 2 + a x di — b-^Ci — a\. Substitute these 
results in 16. On the left-hand side of the equations in 16, since Y is holomorphic in |/| < 1, 
Y x Y~ l and Y t Y~ l can be written as power series in /. And since S is holomorphic in |/| > 1 with 
S(x, t, oo) = J, SxS' 1 and SfS' 1 can be written as power series in 1// with no constant term. 
Thus from the non-negative powers of / in 16, after substituting 18-19, it follows that 

y;y " = (,_„ I) • (20) 

Y.Y-' = (: M - B l + '- i2 ). (21) 
\l 2 -vl-V bl + B y ' 



If X = Y X Y 1 , T — Y t Y 1 then X, T must satisfy the zero-curvature equation 



X t -T x + [X,T] = . 



(22) 



Substituting the forms lpl-lp2 into the zero-curvature equation, required to be true for all /, gives 
the following system of equations: 



v 



b x + b 2 , 



(23) 



B 



\b xx + bh 



X ; 



(24) 
(25) 
(26) 



X 



>t 



(\b xxx + \b 2 x + bb xx + b 2 b x ^ , 
\b xxx + 2b 2 x + bb xx + b 2 b x - V . 



The third equation can be integrated to give V = \b xxx + bb xx + \b 2 x + b 2 b x — S(t), where 5 is an 
arbitrary function of t alone. Using this in the last equation gives 



\b xxx + \b 2 x + 6(t) . 



(27) 



All this can be summarized in the following result: 

Proposition 2.1: Let Z7 (0, 0, Z) be an element of LGL 2 (C), let S^Y be the Birkhoff decom- 



component of l/l in the 1,2-entry of S. Then b(x,t) is a solution, possibly with singularities, of 
the potential KdV equation pkdv for some function 5. 

The reason for the phrase "possibly with singularities" here is because for some values of x and 
t, U(x,t,l) might leave the dense open set of LGZ^C) where Birkhoff decomposition is possible 
(it can be proved that these values are isolated). It is important for the purposes of this paper to 
note that although the above proposition makes no mention of the linear constant-coefficient flow 
linl, the heart of its proof is that the this flow induces, via Birkhoff decomposition, the matrix 
Lax pair lpl-lp2 for the potential KdV equation. Note also that the first equation of the Lax pair 
Ipl gives the usual relation of KdV with the Schrodinger equation. Writing either column of Y 



There are many applications of the above result, of which only one will be discussed here, 
the construction of the standard Backlund transformation for potential KdV. The idea behind 
this Backlund transformation is as follows: Suppose the element U(0, 0, 1) of the loop group gives 
solution U (x, t, I) of the linear system linl with Birkhoff decomposition S* -1 Y and corresponding 
potential KdV solution b(x,t). Let us now try to find the potential KdV solution corresponding 
to the element 






(28) 
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with < 9 < I. The new solution of the linear system linl is 



= rrO o) s " 1( ^ m) ^ m) (i v)- < 29) 

To perform the new Birkhoff decomposition, a certain matrix and its inverse must be inserted as 
follows: 

li-e fo i\ „_-,. . / ol ( a \\_, _ /0 j^' 



S-\x,t,l){ , „ , „ 1 • , „ , „ JF(i,M) i . (30) 



The aim is to chose a and (5 so that this is written in Birkhoff factorized form, i.e. so that 



is holomorphic in |/| > 1 and satisfies S(x,t, oo) = /, and 

Y(x,t,l)=( a 1 )Y{x,t,l)(° ^ (32) 

V ' \l-6 + a(3 (3) V '\l J V; 

is holomorphic in \l\ < 1. Inserting the expansion sexp in scon, the former condition requires 
(3 = b. For the latter condition it is just necessary to check Y does not have a pole at I — 9, and 
this requires a = —Y 2 i(x,t,9)/Yn(x,t,9). Finally, it is necessary to compute the new solution of 
the potential KdV equation, i.e. the component of 1// in the 1,2-entry of S. A brief calculation 
shows this is simply —a. Utilizing the Lax pair lpl-lp2 it is straightforward to determine properties 
of a leading to the following result: 

Proposition 2.2: If b is a solution of the potential KdV equation pkdv and ip satisfies 

ipxx = (9- 2b x )ip , tp t = -\b xx ^ + (9 + b x )ip x (33) 

then b = b + ip x /ip is also a solution of potential KdV, for the same function 5. 

Equations lp comprise the standard scalar Lax pair for the KdV equation. Applying the 
Backlund transformation to the x- independent solution b(t) = J 5(t)dt gives the 1-soliton solutions 

b( x , t) = J 8(t)dt + v^tanh (V9(x + 9t) + C) (34) 

and the singular solutions 

b(x, t) = J 8(t)dt + v^coth (V9(x + 9t) + c) (35) 
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where in both formulae C is a constant. The easiest way to apply the Backlund transforma- 
tion again to this solution is to use the Bianchi permutability theorem that states that the two- 
parameter family of solutions obtained by applying first the Backlund transformation with param- 
eter Oi and then the Backlund transformation with parameter 9 2 is the same as the two-parameter 
family of solutions obtained by applying the two Backlund transformations in the reverse order. 
See mel for a detailed discussion of this. The Bianchi permutability theorem can be used to derive 
an algebraic expression for the solutions obtained by applying two Backlund transformations (see 
dj section 5.4.3): 

Proposition 2.3: If b is a solution of the potential KdV equation pkdv, and bi and b 2 are 
solutions obtained by applying Backlund transformations with parameters 6\ and 9 2 respectively 
to b, then 

B = b+ iFr < 36) 

0i - b 2 

is a solution obtained by applying the two Backlund transformations successively to b, in either 
order. 

Applying this result using a 1-soliton solution for b±, a singular solution for b 2 and 9 2 > Q x gives 
the 2-soliton solution 

J v^itanhai - V^cotho^ 1 a 2 = \[9 2 {x + 9 2 t) + C 2 

= j 5{t)dt + y^itanhai + y^tanha^ (37) 

^tanha^ sech 2 ^! + v^i^tanhcti sech 2 o;2 
\[Q 2 — y/Oi tanhai tanha 2 

This concludes our presentation of the basic theory of the KdV equation and its relation with 
the linear system linl which will be imitated for discrete systems in later sections. 



3 Discretizations I: Lattice KdV 



The aim in this section is to follow the procedures of the last section as closely as possible, but 
replacing the solution U(x,t,l) of linl by the solution U nm {l) of the lattice equation disci, which 
has general solution 



I + h 



r 

1 o 



I + k 



1 
/ 



1 / (1 + hV~l) n + (1 - hV~l) n 

4 v Vi ((i + hViy - (i - hViy 

( (1 + kV~l) m + (1 - kVl) m 
V V~l ((1 + kV~l) m - (1 - kV~l) m ) 



U ,o(l) 

± ((i + hViy - (i - hvir) 
(i + hViy + (i - hViy 

±((i + kV~i) m -(i-kVi) 

(1 + kV~l) m + (1 - kV~l) m 



(38) 



^o,o(/) • 



Suppose U nj7n {l) has a Birkhoff factorization S n ^ n {l)Y rijm {l). Substituting in disci and rearranging 
gives 

1 



V V -1 — <? 



n+l,m 



- 1 n,m+l- f n ,m u n,m+l 



I + h 
I + k 



J 0, 

r 

1 o 



Writing 

1 / On,m b nm \ 
S n , m = I + 7 ( +-.. 
" \ £n,m u>n,m / 

and comparing non-negative powers of / on both sides of intl-int2 gives 

/ 1 - hb n+hm h 

\hl + h(d n+ i !m — a n>m ) 1 + /i6 n>m 

(1 — kb n , rn+ i k 
kl + k(d nirn +i — a„ )m ) 1 + kb nm 



(39) 
(40) 

(41) 



n..m. • 



(42) 
(43) 



There is one further simplification that can be made in these equations, intl (and similarly int2) 
can be written in the form 



c c-i y V~ l — 



I + h 



r 

1 o 



5-1 



(44) 



Taking the determinant gives 

det (5 n , m ^ lim ) det (Y n+hm Y-^ = 1 - hH . 



(45) 



The Birkhoff factorization theorem applies for scalars (elements of LGLiC) too, so from this it 
can be deduced that det (Y n+ i^ m Y~J^j = 1 — h 2 l (and det {Sn,mSn+i,m) = !)■ Applying this to 
int3 (and the corresponding result det (Y n ^ m+1 Y~^j = 1 — k 2 l to int4) gives 



v 

- 1 n+l,m 
Y ri: m+1 



1 — hb n+ i tm h 

hi + &n jm ^n+l,m ^n,rrt^n+l,ra 1 4" hb n ^ m 

1 — kb n , rn+ i k 

kl + fe^m &ri,m+i kb n rn b n m +i 1 + kb n ^ m y 



Y n 



Y 



n,m ■ 



(46) 
(47) 



Up to a rescaling this is precisely Nijhoff ei a/.'s scalar Lax pair for the lattice KdV equation 
Nijhoff. Writing 



1 — hb n+ i >m h 
/i/ -|- b n rn &n+i,m, hb n ^ m b n ^.\^ m 1 -|- hb nr 

A;/ + fo nim b nm +i kb n ^ m b n ^ m j r \ 1 + kb nn 



(48) 
(49) 
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equations dlpl-dlp2 are just 

Yn+l,m L nm Y nrn , Yn,m+1 ^n,mYn,m i (50) 

and for consistency L n ^ m+1 M n ^ m = M„ +lim L n m . (This last equation plays the role of the zero- 
curvature equation in the continuous case.) Substituting the forms found for L n>m ,M njrn in the 
consistency condition gives lattice KdV dp2. Thus the analog of Proposition 2.1 is obtained: 

Proposition 3.1: Let t/o,o(0 be an element of LGL 2 (C), let S~^ m (l)Y n>m (l) be the Birkhoff 
decomposition of 



£Ai,m(0 



/o r 



/o r 



Co,o(<) , (51) 



and let 6„ iJn be (—1) times the component of 1/7 in the 1,2-entry of S n ^ m . Then b n ^ m is a solution, 
possibly with singularities, of the lattice KdV equation dp2. 

In fact there is no reason why £/ ,o(0 should not, in this case, be dependent on h and k. So 
in principle the class of solutions of lattice KdV occuring this way is much larger than the 
corresponding class of solutions of (potential) KdV. 

Let us attempt to find a Backlund transformation and soliton solutions for lattice KdV pro- 
ceeding as in section 2. Making the replacement 



UoM^f^C, oK«(i t) (52) 



gives 



and 



tf»,m(0^f-p(° JW™«(° 'I 9 ) (53) 



V < - V \l-9 + a rhm b n ,m K, m ) \ 1 / 

y,ud - f, """" , ./ t) (55) 

where a n ^ m = — {Y n ^ m )2i{9) / (Y n ^ m )ii{6). The new solution of lattice KdV is simply — a„ im . Using 
dlpl-dlp2 to find properties of a ntm gives the Backlund transformation: 

Proposition 3.2: If b 1hm is a solution of the lattice KdV equation dp2 and i/j n ,m satisfies 

^2, m -2^ B+1>m + ^ B>m = e ^ rhm -( h ^ 2 ' m ~ bn ' m )^ 1 , m (56) 



h 2 ^ ^ 



k h 
10 



+ (b n +l,rn ~ &n,m+l)V'ra,m (57) 



then b ntm = b n+1>m + (Vwi,m - VVi,m) / '(^n,m) is also a solution of lattice KdV. 

The first equation here is a natural discretization of the first equation in lp, and is the discretiza- 
tion of the Schrodinger equation studied in BPPS. The second equation is, however, completely 
unrelated to that in lp. To get 1-soliton and singular solutions the Backlund transformation can 
be applied to the trivial solution 6„ im = 0. This gives solutions of the form 

K m = + ^ + ^ ~ ^ ~ ^ ' kV °) m (A, B constants). (58) 

If A : B is positive and h,k < 1/V9 this gives 1-soliton solutions 

b n ,m = v^tanh (ntanh^Ov 7 ^) +mtanh~ 1 (kV9) + C) (C constant). (59) 
If A : B is negative and h, k < 1/ gensol gives singular solutions 

b n ,m = coth (n tanlT 1 (/i v 7 ^) + m tanh" 1 ( fc v 7 ^) + C) . (60) 

The Bianchi permutability theorem applies equally here in the discrete case, and this can be 
used to give the discrete version of Proposition 2.3: 

Proposition 3.3: If b nym is a solution of the lattice KdV equation dp2, and b^ m and bffl m are 
solutions obtained by applying Backlund transformations with parameters 9\ and 9 2 respectively 
to b n>m , then 

B n ,m = ^>n,rn + ~7J) (2) (^1) 

On,m On,m 

is a solution obtained by applying the two Backlund transformations successively to & n , m , in either 
order. 

Proof. Writing q n ^ m = (ip n +i,m — i>n,m) / (hil> n ,m) , the Backlund transformation can be written as 
b n ,m — > b n+ i tm + q n ,m where q n ^ m satisifies the discrete Riccati equation §s 

(1 - hb n+2 , m + hb n>rn ) + (hO - 

?n+l,m — 7— T (O^J 

1 + hq n>m 

(for the sake of brevity I only look at the first equation in dsll-dsl2). Alternatively, after a little 
algebra, the transformation can be written fe n m — > fe n m where fo n , m , 6„ im are related by 



— # + (&n,m — ^>n+l,m)v>n,m ~ b n +l,m) ■ (63) 



Using the Bianchi permutability theorem and the premises of the theorem, it is known that 
applying the BT with parameter 0\ to b n ^ m gives b^ln, applying the BT with parameter # 2 to b n;m 
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gives b^ m) and applying either the BT with parameter 9 2 to b^ m or the BT with parameter 9 1 to 
bnjn gives the same solution B n , m . This implies 4 equations: 

= 1 + (b n , m - &£> 1>m ) ($1 - b n+1 , m ) , (64) 

= e 2 + (b n , m - bSl^j (b£l - b n+1 , m ) , (65) 

= 2 + (b$> m - B n+1 , m ) (B n>m - b^J , (66) 
= 0i + (&i 2 l - S n+1 , m ) (fl B>m - fc^J . (67) 



1,(1) 

u n+l,m 


- &Sln + 6n+l,m 


bn,m 








l(2) 

w n+l,m 


- + 6n+l,m 


bn,m 








B n +l,m 


R 4- 

D «,m ~r u n+l,m 










B n +\ jjn 


- B + ft (2) 

"ti.m ~ u n+l,m 


_ft(2) 





Adding the first and last of these equations and subtracting the other two gives 

2(^ - 2 ) = (S B>m - 6 n , m )(6gL - &$J + - W,m)(&£i, m - 62 1>m ) (68) 

The general solution of this is clearly 

(S„,m " b n , m )(bW m - bgj = {9, - 9 2 ) + (-l) n F(m) , (69) 

where F is an arbitrary function of m. Using the second equation in dsll-dsl2 is is possible to 
show that F(m) = 0. •. 

All that remains to do in this section is to briefly discuss the nature of soliton solutions of 
lattice KdV, and in particular how they compare to those of continuum KdV. From dsol the speed 
of the soliton with parameter 9 is 

_ htajihr 1 (k\fd) 

(The formal definition of the "speed" is the number c such that the solution depends on m, n only 
through the combination nh + cmk.) Recall that the parameter 9 is limited by the requirements 
hVd,kVd < 1. Thus: 

Proposition 3.4: For h = k the soliton solutions of lattice KdV dp2 all have speed 1. For h < k 
there are solitons with all speeds greater than 1. For h > k there are solitons with all speeds 
between and 1. 

Proof. The result for h — k is obvious. Switching h and k switches c and 1/c so it is just 
necessary to check the result for, say, h < k. As 9 tends to c tends to 1, and as 9 tends to 1/k 2 
(which is less than 1/h 2 ) c tends to oo. So the result will be proved if we can establish that c is 
a monotonic increasing function of 9 for < 9 < 1/k 2 . Writing z = ky/9 and a = h/k < 1, 

atanh -1 z 



tanh 1 az ' 
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(71) 



and it is necessary to check this is a monotonic function of z on < z < 1 for a fixed between 
and 1. Differentiating gives 

T, = M -V)FM(i- AJ ) [ (1 " aV > tanh ~' M " a(1 " - ?) unh " (2) ] ' (72) 

All the terms except the last are evidently positive. The last term can be written g(az) — cxg(z) 
where g(z) = (1 — z 2 ) tanh" 1 (z) . Thus it is necessary to show g(az) > ag(z). But this follows 
immediately from the convexity of g, which is trivial as 

d 2 a 1z 

— | = -2 tanh- 1 z < for < z < 1 . • (73) 

dz 2 1 — z 2 

Proposition 3.4 does a lot to clarify the relationship of lattice KdV dp2 with its standard 
continuum limit b t = b x on the one hand, and potential KdV pkdv on the other. The linear 
equation b t = b x admits solitons of speed 1, but, since it is linear, the solitons can be of arbitrary 
amplitude. The indirect method of discretization used has given rise to a nonlinear discretization, 
except when k = h (when dp2 can be written as a product of linear factors). The family of speed 
1 solitons with arbitrary amplitude is perturbed, after discretization, into a family of solitons 
with a nontrivial speed-amplitude relation. For small h, k the low amplitude solitons (those with 
\f~9~ <C l/h,l/k) must have speed close to 1, and indeed this is the case. For larger amplitudes 
the speeds can change substantially, giving a range of speeds ranging from 1 to either or oo. 
Since there now are solitons of different speeds, and the necessary algebraic structure has been 
preserved, the phenomena associated with KdV will emerge, in particular elastic soliton scattering. 
Thus from a phenomeno logical viewpoint, lattice KdV is closer to potential KdV than the linear 
equation bt = b x . There are, however, several fundamental differences: First, the range of soliton 
speeds in lattice KdV is limited to speeds either less than or greater than 1. Second, there are 
many solutions of lattice KdV that do not have natural continuum limits; for example, solutions 
gensol in the case where 9 exceeds 1/h or 1/k (or both). 



4 Discretizations II: The simplest natural discretization 

This section is devoted to the simple discretization disc2 of linl, which, as explained in the 
introduction, should give an integrable lattice equation which has potential KdV as its standard 
continuum limit. The general solution of disc2 is given by huge on replacing k with kl. 

Once again suppose U n , m {l) has a Birkhoff factorization S , ~j n (/)V„ irrt (/), and substitute in disc2 
to get 

Sn' m , (74) 



1 n+l,™ 1 n ,m 



= s„ 



n+l,m 



I + h 



1 

1 
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Y, 



n.m+l-f n ,m 



y-1 



m+1 



I + kl 



o r 

Z 



5-1 



Writing 



1 / Q"n,m bn,m 
Jn,m = I + T I 



I V Cn in G?^, 



1 / Q j n,m 



+ ... , 



(75) 



(76) 



n, m •- u n,m , 



and employing the relations det (Vn+i^Y^) = 1 — /i 2 /, det (V^m+iY^) = 1 — /c 2 / 3 (obtained by 
left-multiplying tml and tm2 respectively by S ntTn S~l ltm and S n , m S~^ m+1 , taking the determinant 
and factorizing) gives 



where 



V - T V V — M V 

1 n+l,m ± -<n,m 1 n,m i 1 n,m+l lvl n,m 1 n,m i 



1 — hb n+ i )m h \ 

hi ~\~ b n rn ^n+i,m hb n ^ m b n -\-\^ m 1 -|- hb n rn J 

( ^>n,m-\-\ -\- Yu nm -\- A^^ j3n,mbn,m+l 



M, 



n,m 



= k 



A + Pn,rn \ 



/^n rrt ^n,m(bn,m ^n,m+l) "I - ^- l n,m(bn,m ^n,m+l) ^n,m "I - Pn,mbn,m 



(77) 



(78) 



.(79) 



The matrix M depends on three lattice fields f3, A, E in addition to the basic lattice field 6, but 
£ is determined via the relation 



y 



+ A 2 + 8 3 . 

' n,m 1 r"n,m 



k 2 



(80) 



Substituting these ansatze into the consistency equation cons gives the following 3 equations for 
the 3 fundamental fields b, (3, A: 

bn+l,m bn+l,m+l ^n,m ^n,m+l 



+ {b n +l,m+l ~ b nm )(b njm+ i — & n +l,m) , (81) 



A n+ i ;Tn + A n>m — f — - : j f — 1 + — (fc n +i, m +i + & n +i,m ~~ b rhm+ i — fe n , m ) j , (82) 



h J \ k J 

\A ^ 2 (^n+l,m + /^n+l,m) ~~ ^1 + /^(A 2 ^ + /? 3 m ) 



i/i/c 2 



(83) 



Note the equations involve 6 at 4 points (& n , m , & n +i,m, b n>m+ i, & n +i,m+i) but /3 and A at only 2 
(/3n,m,/3n+i,m,A nm ,A n+1)m ). The system fkdvl-fkdv3 will be given the title /u/Z lattice KdV; as 
will shortly be shown, unlike standard lattice KdV, full lattice KdV displays, for certain choices 
of h and fc, solitons with the full range of speeds. Full lattice KdV also has, as expected, potential 
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KdV as a standard continuum limit: Replacing b n , m by b(x,t), b n+ i >m by b(x + h,t), & n ,m+i by 
b(x,t + k), b n+ i tm+ i by b(x + h,t + k), and similarly for (3 and A, and then taking the limit 
ft,, &; — > 0, the equations fkdvl-fkdv3 become 

2(3 = 2b X} 2A = -(3 X , 2(3 x b t = (A 2 + (3 3 ) x . (84) 



Eliminating (3 and A from these yields potential KdV b t = \b xxx + |6 2 



2 x 



There are analogs for full lattice KdV of all the results of the previous sections: 

Proposition 4.1: Let U 0t0 (l) be an element of LGL 2 (C), let S~^ m (l)Y n>m (l) be the Birkhoff 
decomposition of 



I + h 



o r 
/ o 



n 



I + kl 



1 
/ 



Uo,o(l) , (85) 



and let b n ^ m be (—1) times the component of l/l in the 1,2-entry of S n ^ m . Then fe n m is a solution, 
possibly with singularities, of the full lattice KdV system fkdvl-fkdv3. 

By "6 n ,m is a solution of full lattice KdV," I mean that there exist fields (3, A for which equations 
fkdvl-fkdv3 hold. In practice, once b is known, the easiest way to determine (3, A will be directly 
from equations fkdvl and fkdv2. In the previous proposition the other fields can actually be 
determined from S if this is known in full: If the expansion of S in powers of l/l is as in S, then 



ftn,m ^n,m+l d nm b nrn b nm +\ (86) 

1 

2 



(bn,m+l + ){a 

n,m + a n,m+l ^n,rrt&n,m+l) ) (87) 



k 2 

+ (^n,m+l — &n,m)(o n ,m + Qn,m+1 ~~ ^n,m&ra,rrt+l) ) (88) 

The Backlund transformation takes the following form: 

Proposition 4.2: If 6„ im , /3 n , m , A n>m is a solution of the full lattice KdV equation fkdvl-fkdv3 
and satisfies 

^2 = ^n,m - I ^ j V>n+l, m (89) 

^ = l^ + Pn, m j ^ (yuj 

+ ^(^n+l,m — &n,m+l)($ + Pn,m) + A Ttjm + S„ im — — ^ ^Vi,™ 

then 6°^™ = b n+ i iTn + (ip n +i,m — VVi,m) I (hip n ,m) is also a solution of full lattice KdV. The fields /3, A 
are replaced by /3 new , A new respectively, which are given by the following algebraic equations: 

Pn°Z+Pn,m = + ( 6 n , m - 6^ +1 ) (6JJ^ - &n,m+l ) > ( 91 ) 
^n,m + An,m = — (/3 nm — (3 n ^ (b n rn+1 + 6„ m — & n , m +l — & n ,m) • (92) 
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The Backlund transformation for continuum KdV (Proposition 2.2) is recovered in the limit h,k — > 
0, since /3 njTn — > b x , A„ iTO — > — \b xx , S„ im — | — > 0. Note the difference between the second equa- 
tion in fsll-fsl2 and the discrete evolution proposed in BPPS. The solutions obtained using the 
Backlund transformation on the vacuum solution b ntin = (3 ntm = A„ iJn = are given by gensol with 
k replaced by k6. In particular, writing t(n, m) in place of tanh (n 
it is straightforward to verify that the soliton solution is given by 

b n ,m = V0t(n,m) , (93) 
t(n,m + 1) -t(n,m) 

t(n,m + l) 2 -t(n,m) 2 

1 (t(n,m + 1) -t(n,m)f 
"' m = 2k ' ( > 

Before exploring the phenomenology of these solitons, note that since the proof of Proposition 
3.3 is based almost entirely on the first equation of the scalar Lax pair dsll-dsl2, it is no surprise 
that it goes through verbatim to full lattice KdV, i.e. 

Proposition 4.3: If b rhm is a solution of the full lattice KdV system fkdvl-fkdv3, and b$ m and 
b^ m are solutions obtained by applying Backlund transformations with parameters 9\ and 6 2 
respectively to 6 n ,m, then 

B n ,rn — b n ,rn + 7j2)~ (^) 

On,m On,m 

is a solution obtained by applying the two Backlund transformations successively to b n>m , in either 
order. 



Since the formulae algl-alg2 for applying the Backlund transformation to the fields f3, A are 
already pure algebraic there is no need to consider them in proposition 4.3. 

It just remains to investigate the speed-amplitude relation of the soliton solutions. The soliton 
speed is 

_ htanh'^keVd) 
°~ ktanh^ihVO) ' 

where the range of the parameter 6 is limited by the requirements k6y/9,hy/6 < 1. Writing 
a = hk~ 1/3 and z = V9k 1/3 gives 

, 9/ ott tanh _1 (z 3 ) , . 

c = k 1 ; -i, r • (") 

tanh [az) 

See figure 1. For a < 1 <^ h 3 < k the speed is a monotonic increasing function of z (or 9), going 
from as z — > to oo as z — > 1. Thus for this range of parameters the soliton content exactly 
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Figure 1: The function a tanh _1 (2 3 )/tanh _1 (a^) for a = 0.1,0.8, 1, 1.02, 1.2, 1.5,2.0 (from top to 
bottom) 

mirrors that of continuum KdV. For a = 1 <^ h 3 = k the speed is a monotonic increasing function 
of z (or 8), going from as z — > to 1 as z — > 1. For a > 1 <^> /i 3 > there is an interesting effect 
that c increases from 0, reaches a maximum value, and then decreases again to as z approaches 
I /a. Thus for these choices of h, k there is a limited set of speeds, but for all but the fastest there 
are solitons of 2 different amplitudes; furthermore these can be superposed to give other types 
of soliton solution. Note that if our interest in discretizations of KdV were for the purposes of 
numerical simulation, we would presumably want both h and k small and of the same order of 
magnitude, and thus be in the h 3 < k regime, where the soliton phenomenology is correct. 

5 Discretizations III: A second-order discretization 

In this section our method is applied to the discretization sol-so2 of lin2 with N — 1. The resulting 
system is of limited intrinsic interest, the main point here is to illustrate that our methods can 
in principle be extended to give a whole range of integrable discretizations of equations in the 
KdV hierarchy. One interesting point that emerges is the form of the related discretization of the 
Schrodinger equation. 

Following the usual procedure, assuming U n ,m(l) has a Birkhoff decomposition S~^ m {l)Y ntm (l) , 
writing 




(100) 
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etc., gives the system 



V —TV V — M V 

1 n+l,m J - J n,m 1 n,rn i 1 n,m+l lvl n,m 1 n,m i 



(101) 



where 

L 7) ™ 



1 + \h 2 l — hb n+ i tJn + \h 2 ( y a n+ i, m — a„ jm ) h + \h 2 {b n , rn — 6 n +i jTn ) 

/iZ + b n ^ m — b n+ \. m — hA 1 + \h 2 l + hb ntTn + \h 2 (a n ^ m — a n +i 5 

(1 — kb nyfn+ i k \ 

kl ~\~ b n ^ m &n,m+i kb n ^ m b n ^ m +\ 1 -|- kb n ^ m ) 



and 



2 (^ri,m ^n+l,m) 2 ^'( G! 'n.+l,m. ^n,m)(^n+l,m ~l~ b nrn ) 4~ ^/l (d n ,m a n+l,m) 



(102) 
(103) 

(104) 



1 + \h{\) n)m — & n +l,m) 

The consistency condition L n ^ m+ iM n ^ m = Mn+i^L^m unravels to two equations for the fields a, 6. 
Introducing the combinations 



J n,m 



^n,m 



2 (^n+l,m+l a n,m+l 4" On+l,m Q"n,m) i 
1 



the equations can be written 
1 ^6. 



■ J n,m 

= 

+ 



/i 2 

+ T 



(^n+l,m+l ^rt,m+l ^n+l,m "I - ^n,m j r h h 

j£ r On+l,m+l°n+l,m ~~ °n,m+lO n 

2 _J ^ £^ (^n+l,m+l + ^n,m+l ~~ &n+l,m ~~ ^n,m) 

^ (^n+l,m+l ~~ frn,m+l + ^n+l,m — ^n,m) 

(^n, m(^n+l, m+l^ra, m ^rt+1, m^n, i 



.m ? 



'/i 2 A 



(105) 
(106) 

(107) 
(108) 



i+l) — ^n,m ~~ ^n,m) ~~ 



+ ^ (&n,m+l + &n+l,m — m ~ &n+l,m+l + 4&ra+l,m+l&n,m _ 46 n+ i 5 m & nj m+1 ) 



4A; 
h 2 

+ -T- 



/i 2 /i — h 2 h 4- /i 2 h — h 2 h 

u n+l,m+l u n+l,m u n,m u n,m+l ~r u n ,m+l u n,m ^ n +i )m +i" 



+ 4 



8 ^ V - & n+l,m & n+l,m+l + &n+l,nA,m+l ~ ^,m+l fo n+l 

1 ^ ^ . . A .. _L% A h . . .. h 

( 



n, m 



, m^n+l,m+l 3& ra ,m^n,m+l ^rt+1, m+1 ^n, m \ 
2 \ — 2b rh m b n+ i : m — & n +l, m m+1 — 26 n+ i ) m+ i 6 Ttj m+ i / 



V-6 2 

,2 



/l 2 

+ y 



1—1 . „ j I ' 7 ' '"3 ' ' 

3&n+l,m+l^n+l,m^n, m+1 + 36 ri +i j m & n +l, m+l^n, m 
— 36n, m+l^n, m^n+1, m+1 36 n _|_i ) mb-n, mbn, m+1 

jAt+l,m — ^n+l,m^"+l,m+l + ^n,m^n,m+l + ^n,m+l^n,» 

^n,m+l^n,m + ^n+1, m^n+1, m+1 + 26 n +l, m^n+1, m+l^n, m&n, m+1 ^ 

~ ^n+l,m^n+l,m+l^ri,m+l ~ &n+l, m+l^n, m V m+1 
V — &7j, m+l^n, m.bn+1, m ~ ^n+1, m^n+l, m +l^n, m / 
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Since the field a only appears in the equations through the combinations A and E, which only 
depend on a through differences, solutions of this system are only defined up to addition of a 
constant to a. The analog of propositions 2.1, 3.1 and 4.1 is 

Proposition 5.1: Let U 0j0 (l) be an element of LGL 2 (C), let S~^ m {l)Y n>m {l) be the Birkhoff 
decomposition of 

(l + OL h \ n /l k\ m 

u nA0 = , , 2 , , hH , , J u 0fi (l) , (109) 



hi 



1 + — 



and let a„ im and 6 nm be, respectively, the 1,1-entry and (—1) times the 1,2-entry in the 1/7 
component of S njTn . Then a„ jm , 6„ i?n is a solution, possibly with singularities, of the system 
mdkdvl-mdkdv2 . 

The system mdkdvl-mdkdv2, despite its algebraic complexity, is an integrable discretization 
of the equation b t = b x in every sense that lattice KdV is. The soliton solutions are given as 
follows: 

Proposition 5.2: The system mdkdvl-mdkdv2 has soliton solutions 

bn,m 



V0 tanh (n tanh -1 ^ - ^^ 2g j + m tanh~ 1 (kV9) + C^j , (110) 

(111) 



= constant, 



with speed 



i ( hVe 



/ctanh 



(112) 



In greater generality, it can be shown that if instead of equation sol a pth order approximation 



U, 



n+l,m 



u„ 



(113) 



is used, then the speed of the soliton solution is 



A; tanh" 



-i / s 



,(hVe) \ 



(114) 



c p (hVe) J 

where c p (x) and s p (x) are, respectively, the order p truncations of the Taylor series for cosh(a;) 
and sinh(:r) (ignoring terms of order x p+1 and higher). It is straightforward to verify that for 
small x 

' - ,= S + 0(^) xodd (U5) 



X \ ^ X J J 

Thus for small h the dependence of the soliton speed on h becomes weaker as p increases. Likewise 
the order of accuracy in k can be increased. (The distinction between the even and odd cases in 
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tanap, that for odd p there is a "free" extra order of magnitude accuracy, means that mdkdvl- 
mdkdv2, for which p = 2, is actually no more accurate in this regard than standard lattice KdV, 
with p — 1. The equations obtained from p = 3 can be written down, but due to their length I 
have restricted the discussion to the p = 2 case.) 

Returning to the formula hewspeed, note that if v — ky/Q, w = hy/0, then 

- tanh -1 v 



i tanh -1 



(116) 



The function in the numerator increases monotonically from 1 to oo as v goes from to 1. The 
function in the denominator decreases monotonically from 1 to as v goes from to oo. Thus 
for the current discretization c can only take values greater than 1. 

The soliton solutions just presented can be found using the Backlund transformation, which 
is obtained as in previous sections: 

Proposition 5.3: If 6 n ,m, a n,m is a solution of mdkdvl-mdkdv2 and ip n ,m satisfies 

VVi+2,m - 2VVi+l,m + VVi,m _ q , I , h ■ \ 



+ 



n+l,m 



(b n +l,m ~ b ntm )ll> n+ 2,m ~ 3(6 n+ 2 n, 

2h 

h 2 6 2 I h(b n+ i m — b n+ 2,m) \ i a n+2,m ~ ^ a n+l,m + Q n ,m , 

J- ^ + 2 J Vn ' m + 2 ^ 

(&n+2 jn + )b n +i 
+ ^ n+i,m 

/i / Q"n,mbn+l,m ~\~ Q"n+1 ,m,bn+2 ,m Q"n+2,mbn,m \ 
+ 7 _ , , 7 1 V'n+l.m 

^ \ ^n,m"n+2,m Q>n+l,m"n,m ^n+2,m"n+l,m / 



h 



+ (&n+l 



1pn,m+l j>n,m \ — f (0 + a n +l ;m — a n ,m + ^n,m+l(^rt,m ~ ^ra+l.m)) VV, 

1 + 2^(^n,m ~~ ^n+l,m) 



(118) 



then 



7 I V'n + l.m-l/'n.m , lu2( n _ _ _ /1\ 

°ra,m — -, . lr/7 7 \ l iiy J 

<:Z = b n , m b™ - a n , m + constant (120) 
is also a solution of mdkdvl-mdkdv2. 

All formulae in the previous proposition have been written in a manner that hopefully makes it 
clear in what sense they are modifications of the corresponding formulae in proposition 3.2. The 
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surprising feature of the discretization of the Schrodinger equation in proposition 3.2, equation 
dsll, is that in it the parameter 9 multiplies ip n ,m, not Vwi,m, which would seem more natural. 
The new discretization just presented, equation mdkdvl, has 9 multiplying ip n +i tm . But the cost 
of this is the introduction of many new terms, including a term proportional to 9 2 , multiplying 
V\i,m- It can be checked that the new discretization mdkdvl is a second order approximation to 
the Schrodinger equation, while dsll is only first order. This is the justification for the title of 
this section. 

6 Concluding Remarks 

In this paper I have presented a systematic approach towards integrable discretizations, based 
on the loop group approach to integrable systems. Three integrable discretizations have been 
examined in detail, one known, the lattice KdV system of Nijhoff et al, and two new, one of 
I have called full lattice KdV, as it would seem to be the first discrete integrable system with 
(potential) KdV as its standard continuum limit. For each integrable discretization a Backlund 
transformation has been given and soliton solutions have been analyzed. Unlike the lattice KdV 
system of Nijhoff et al., which only displays solitons with speeds below, above or equal to 1, full 
lattice KdV has the full range of soliton speeds (for suitable choices of h and k). 

Full lattice KdV would seem to merit further attention. Our plans include conducting numer- 
ical studies, and to try to work out a suitable inverse scattering formalism. Another issue that 
has not been touched upon in this paper is the subject of tau functions for discretizations. The 
linear flows on a loop group that underlie KdV can be extended to the central extension of the 
group, and one would expect the same to be true for the discretizations looked at in this paper. 

The formalism developed here can also be extended to look at integrable discretizations of 
KdV on non-rectangular lattices, see hex. 
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